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Development of an aneutronic Fusion Reactor 
Based upon the n B (p,a)aa reaction 


Tri-Alpha Energy, Inc. - Foothills Ranch,Ca. 

Controlled fusion in a field reversed 
configuration and direct energy conversion. 

(US patent: 7459654) 

The U B (p,a)aa reaction at low energies is dominated by the 2'T=1 
resonance at 675 keV which has a width of 300 keV (T p = 150 keV, 
r al = 150 keV). 

The Q-value of the reaction is 8.58 MeV (E in(cm) x 14). 
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The n B (p,a)aa reaction 


When 675 keV protons strike n B, a resonance in 12 C is 
formed at 16.6 MeV having a width of 300 keV. 

This resonance decays by emitting a primary a-particle, 
leading to the first excited state of 8 Be which is 3 MeV 
above the ground state. This state decays into two 
secondary a-particles. a(a 1 ) = 600 mb! 

The state-of-the-art understanding of this reaction came 
from a 1987 study which "proved" that the reaction 
yielded one high energy primary a-particle, and two low 
energy secondary ones. 
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Tri-Alpha Energy Corporation 


We have proven that this is incorrect! 

But first— a few words on the concepts which 
underlie the aneutronic fusion reactor presently 
being developed at T ri-Alpha, Inc. 
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Magnetic Field map of the Field Reversed Configuration 
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FIG , 2 
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Annular layer of plasma rotates around the null surface (86). 
Collisions between injected p's and n B create "high energy” 
alpha particles. Relative energy of protons to n B is about 600 

keV. 
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The direct energy conversion system utilizes an Inverse Cyclotron 
Converter (ICC). D's are replaced by four semi-cylindrical 
electrodes. Energy is removed from the alphas as they spiral past 
the electrodes connected to a resonant circuit. 
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The design of a 100 MW reactor is underway. 
Test " shots " to demonstrate plasma 
confinement are in progress. 




Brief History of the n B (p,a)aa reaction 

The history of this reaction is almost as long as the history of nuclear 
physics itself. 

Lord Rutherford studied the reaction over 75 years ago at -200 keV at 
the Cavendish Laboratory, measuring the ranges of particles coming 
from the target. In his 1933 Proceedings of the Royal Society paper 
he states: 

"we might anticipate that the most probable mode of disintegration 
would be for the three a-particies to escape symmetrically with 
equal energies". 

(Pmc. Roy. Soc. Lond.A 141, 259 (1933).) 
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Three years later (1936), Dee and Gilbert, also of the Cavendish, 
published the results of their expansion chamber studies of the 

n B (p,a)aa reaction. Their 300 keV studies led them to conclude: 

"the common mode of disintegration is into two [alpha] particles 
which proceed at angles of 150° to 180° relatively to one another, 
the third particle receiving very little energy" 

The track images shown in their paper strongly support this 
interpretation. 

(Proc. Royal Soc. Lond.A 154, 279 (1936).) 
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This photo shows the typical mode of disintegration of boron into three a- 
particles. A and B are emitted in nearly opposite directions while the third a- 
particle C receives very little energy and barely emerges beyond the beam of 

scattered protons. 

Downloaded from rspa.royalsodetypublishing.org on October 3, 2012 

Dee and Gilbert Proc. Roy. Soc. A , vol. 154, Plate 5 
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Brief History 

Fast forward to 1987. Of course many experimental and 
theoretical works were performed in the interim. But the 
work of Becker, Rolfs and Trautvetter (Z. Phys. A327, 341 
(1987)) established the modem view of this reaction: 

“The reaction mechanisms have been studied via 
kinematically complete coincidence measurements 
showing that the reaction proceeds predominantly by 
a sequential decay via 8 Be” 
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a-YIELD (REL. UNITS) 


Primary and secondary distributions from Becker et al. The 
two are about equal, leading to the factor of " 2 " in 

determining the cross section. 





The u B(p,a)aa reaction 


"We want to know the energy and 
location of every outgoing alpha 
particle " Tri-A Ipha E nergy Inc. 
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Results of measuring the a-yields as a function of 
incident proton and outgoing a-energies. 


Previous results of Becker etal. indicated that only two alphas were 
observed in finite detectors. We find 3. NACRE assumed 2 in its 
Astrophysical S-factor compilation for this reaction. ( This needs 
correcting) 

As will be shown later, the numberof_al£has is a function ofboth the 
Qroton energy and, the alpha energy. 

To avoid this issue, we plot Counts instead of cross section, normalized 
to the integrated luminosity and solid angle. Cross sections can be 
obtained by dividing by the expected number of alpha particles, but 

the absolute number of a's is what is important to Tri-Alpha. 
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Counts/LA£2 [mb/sr] 


Results showing the dominant 0.675 MeV resonance 
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Two-step Reaction Simulation 


Use experimental a-a phase shifts to generate the resonance 
line shape of 8 Be (2+). A pure single level assumption allows 
correction for potential scattering using hard-sphere phase shifts. 

6 2 = - (|>2 + tan' 1 ( r/2 ) 

E r + Ar _ E 

The second term is what we are after. 

(|> 2 is the "hard-sphere" phase shift... easily calculated. 

This avoids having to compute the energy dependence of r and 
A. It all comes from the experimentally determined values of the 
elastic scattering nuclear phase shifts § 2 . 
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Two-step Reaction Simulation 


Weight this with the probability that the a escapes from the 2' 
resonance in the 12 C nucleus (the penetrability for 1=1 (or 3) alpha 
particles) to get the probability distribution for populating the 8 Be* 
first excited state which subsequently decays into two secondary 
a particles having equal energies in the 8 Be* frame. 

Randomly distribute the primary and secondary a's in the center- 
of-mass frame while preserving proper kinematics and boosts 
bac k to the lab frame. 
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Angular distributions 

We need to know the form of the angular 
distribution of both the primary alphas, and 
of the secondary ones wrt. to the axis 
defined by the primary alphas. 

These can be calculated using the Z bar 
formalism of Blatt and Biedenharn. 
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For the 2‘ state at 0.675 MeV 

For the 2- (un-natural parity) state at 0.675 MeV, we 
have 1=0 and l'=l (according to Becker etal.). 
This immediately leads to an isotropic angular 
distribution for outgoing a's, as observed 
experimentally. 

The distribution of the secondary a's wrt. the 
internal primary axis is found to be a(0) = 1 - P 2 
(cos0), which is proportional to sin 2 0. 

This is when the fun begins! 
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F ailure of the two-step model at the 2- resonance at 
0.675 M eV (r= 300 keV). 
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What went wrong? 

When the 2' resonance in 12 C decays to the first excited 2 
state of 8 Be, the primary a-particle can have orbital 
angular momentum ofl or 3 units. 


The 1987 work assumed the value ofl. 1/1/e tried 3! 
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To change to 1=3 outgoing primary alpha particles 
requires two changes in our simulation 

1. Replace the 1=1 penetrabilites by 1=3 ones when 
generating the line shape of the 8 Be 2 + state. 

2. Use the angular distribution for the secondary alpha 
particles when the primary alphas have 1=3. 

This is obtained by considering 
(5=2- + 1=3 -> 2 + -> S'=0 + 1'=2 
From Z bar (3232;2L) Z bar (2222;0L) P L (cos0) 

We obtain: o(0) = C [1.0 + 2/7 P 2 (cos0) - 9/7 P 4 (cos0) ] 
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Angular Distribution [a.u.] 


Which looks like 
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The effect of the 1=3 penetrability on the resonance line 

shape: a much narrower peak! 

(down from 1.5 to ~0.7 MeVN!) 
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Success of the two step model with 1=3 
primary alpha particles 



Light Nuclei from First Principles- 
INT-2012 


And the angular distribution remains isotropic as observed. 


Ep= 0.65 MeV 
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What can additional experiments 

tell us? 

We performed a coincidence experiment using position 
sensitive solid state detectors to increase the solid angle 
and provide data over an extended angular range. 

Two detectors which each subtended an angular range of 
+/- 8° were employed. Measured angular resolution was 
0.2 degrees. 

The detectors were placed symmetrically on the left and 
right sides of the beam. An opening angle range of 100° 
- 180° was scanned. 
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E„ [MeV] EJMeV] 


Coincidence spectra for the 0.675 MeV resonance (top) and the (3 ) 2.64 MeV 
resonance (bottom) indicating two high energy alpha particles at 0.675 MeV. The 

lab a-a opening angle is 150° 
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Testina the model 


The simulation was asked to predict the 
coincidence counting rate as a function of 
opening angle. Coincidence data were taken 
for opening angles ranging from 100°-to-180°. 


INT- 2012 


Counts 


Both a-particles were required to have energies greater than 3 M eV. 
Excellent agreement is observed. 
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TUNL Data 


1=3 Simulation 
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Conclusions 


l/l/e have found that a two-step sequential model with 1=3 
primary a-particles is required to describe the data at 
the 0.675 MeV 2' state. 

The 1=3 assumption predicts the existence of two high- 
energy a-particles at an opening angle centered at 
155° -as originally proposed by Dee and Gilbert in 

1936! Our coincidence measurements confirm this 

result. 
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Huge impact on the reactor design! 


Our discovery of TWO high energy a-particles is 
having a huge impact on the reactor design. 
They are much easier to extract and convert 
more efficiently into electricity. 

Tri-Alpha Energy is thrilled! 

(our results are published in Physics Letters B 
696 (2011) 26-29 and The J ournal of Fusion 
Energy (Springerlink.com), M.C. Spraker etal., 
2011 ) 
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Questions remain... . 

A challenge to theorists— 

Why does the 2 - resonance at 16.6 MeV decay via 1=3? 

This resonance is assigned T=l, based upon its em decay 
properties. Also, the other members of the isospin triplet 
are observed in 12 B and 12 N (second excited states). 

How does this T=1 state decay into 3 alphas? Obviously 
there is strong isospin mixing. 
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Astrophysics 


Triplc-a Process: 



M. Hjarth- Jensen. Physics 4. 38 (2011) 
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Hoyle State 





o + 
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Hoyle State 



o+ 





Fynbo L: Freer. Physics 4. 94 (2011) 


Light Nuclei from First Principles- 
INT-2012 




More recently, Epelbaum etal. (Lattice EFT) find an obtuse triangular 

configuration. 



FIG. 3: This shows Initial state F , a wavefunction consisting 
of three alpha dusters formed by Gaussian packets centered, 
on the vertices of a- bent-arm or obtuse triangular configura- 
tion. There are a total of 24 equivalent orientations of this 
configuration. 
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12 C 
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0 7-ray beam will not excite 
G + state and will excite 3 “ 
state with very small 
probability. 
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"'Be 4- a 



0 7-ray beam will not excite 
0 + state and will excite 3 “ 
state with very small 
probability. 

9 Measuring a-particle 
angular distribution gives 
J~ of state. 
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Measuring 12 C(7, a) 8 Be requires: 
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12 C( 7 ,a) 8 Be 


Measuring 12 C(7*a) 8 Be requires: 
o Intense, monoenergetic 7-ray beam 
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12 C( 7 ,a) 8 Be 


Measuring 12 C{7 T a) 8 Be requires: 

o Intense, monoenergetic 7-ray beam 

9 Detector capable of measuring angular distributions of 
recoiling a -particles with little or no background 
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HlyS 


•High Intensity y-ray Source 

(HlyS) 


-Located at the Duke Free Electron 
Laser Laboratory 
Part of the Triangle Universities 
Nuclear Laboratory (TUNL) 

-Intra-cavity Compton Back 
Scattering of FEL photons by 
electrons circulating in the Duke 
Storage Ring 
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RF Cavity 


Provides circularly and linearly polarized, nearly monoenergetic y-rays from 2 to 100 MeV 

Utilizes Compton backscattering to generate y-rays 


Optical Klystron 

FEL 


Booster Injector 


Mirror 


LINAC 


HlyS y-ray beam generation 
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Two Bunch Mode 


RF Cavity Linac 



Created by Brent Perdue, 2005 
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Some typical beam intensities 


E jMeV) 
1 - 2 
8 - 16 
20- 45 
50 - 100 
^160 


Beam on target (AE/E =3%) 

2 x 10 7 y/s 

2 x 10 8 (total flux of 2 x 10 9 ) 

8 x 10 6 

4 x 10 6 (will increase by x3-5 in 2012) 

1.2 xlO 7 (by 2015) 
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Show HlyS animation. 
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hi 7 s O-TPC 


mirror optics chain 
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hi 7 s o-tpc 


Drift Chamber 


second avalanche grid 
first avalanche grid 
anode grid 


cathode grid 
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hi 7 s O-TPC 


Optics Chain 


electrostatic demagnifier 


micro-channel plate 



"V 

auxiliary lens assembly 


- main lens assembly 
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hi 7 s o-tpc 


mirror 


optics chain 



lOfflcfl 


- photomultiplier tubes 


- drift chamber 
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[cm| 


Background Rejection 



Cosmic 
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Event Identification 


12 C( 7 , a)*Be 16 0( 7 ,q) 12 C 



-10 -5 0 5 10 -10 -5 ft 5 ID 

[cm | [cm | 
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Event Identification 


12 C(7,Q) 8 Be 16 0(7,a) 12 C 
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Photo-electrons 


Time Projection 


Tilted Track 



Earl i est a rr i v a I Latest a riva I 





Time 


Flat Track 



Drift Volume 


Drift Volume 


Time 
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Time Projection 


12 C(7,a) 8 Be 


16 0( 7 ,cv) 12 C 
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12 C(7,«)*Be 


16 0( 7 ,a) 12 C 



LUX ckimifel [lOni^r channel] 
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Time Projection 
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Angular Distributions 


a 0 was calculated for each event from the track image and 
from the time projection. 
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Angular Distributions 


a 0 was calculated for each event from the track image and 
from the time projection. 

9 Angular distributions were fit in terms of \E1\, |1?2| ? 
W{$) = | sill 3 0 

x (3|£’1| J H- 25| J E'2| 2 cos 2 ^ 4- 10 a/3||151 | J E 2 |cos 0 i 2 Cos^ i 
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l 





Angular Distributions 


a 0 was calculated for each event from the track image and 
from the time projection. 

a Angular distributions were fit in terms of |i?l|, |i?2| ? & 12 : 


lh(d) = — sin" O 

Jsal 

x (3|£1| 2 + 25|£2| 2 cos 2 0 4- 10\/3|£1||£2| COS <j>l2 cos 0 


+ Since angular information is available for each event, 
Unbinncd Maximum Likelihood fits were used. 


J C(|£l| ; |£2|> 12 )=p[^ i ) 

i= 1 
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Angular Distributions 
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Cross Section 


_Q 

4 


b 



E 7 (MeV) 
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Cross Section 



Eres r Q (res) r 7o (res) B(£2 : 2+ -4 0, ) 

(McV) (kcV) (meV) (e 2 fm 4 ) 

10 . 03 ( 11 ) 800 ( 130 ) 60 ( 10 ) 0 . 73 ( 13 ) 
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Phase 


E1-E2 phase difference: 

012 = 5% — S\ + arctan ( 17 / 2 ) 
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Phase 


EI-E 2 phase difference: 


<j> 12 = ^2 — + axe tan (77 / 2 ) 


Nuclear phase shifts: 

Si = arctan 
Hard sphere scattering phase shift: 




arctan 
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Phase 
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N 2 (ct v ) (cm 6 mo I 


Triple-a Reaction Rates 



Temperature (ID 9 K) 
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N 2 {cr x) ) ( c m 6 mo I 


Triple-a Reaction Rates 



Temperature (10 9 K) 
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Carla Froehlich (NCSU)— Simulation of ejecta from supernova explosions 
based upon nucleosynthesis via the vp process. Explains the large abundance 
of Sr observed in metal poor stars. Our results (vs J INA) impact the production 

of elements having A>90. 

Ref. Arcones, Froehlich and Martinez-Pinedo, ApJ 750:18 (2012) 
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Summary 


r= 


2 2 


1 

Eras 

r a (rcs) 

r 70 (res) 

B(E2 : 2+ -> Of) 


(McV) 

(kcV) 

(mcV) 

(e 2 fm 4 ) 


10.03(11) 

800(130) 

60(10) 

0.73(13) 



o 2t in lL> C has been directly observed through the 
12 C + 7 — > 3o reaction. 


Light Nuclei from First Principles- 
INT-2012 




Summary 


f 

2 s 


r 

E a (res) 

r 7o( reS ) 

B(E2 : 2} -» Of) 

(MeV) 

(kcV) 

(mcV) 

(e 2 fm 4 ) 

HI 7 S 10.03(11) 
EFT 9.65 

800(130) 

60(10) 

0.73(13) 

2 ( 1 ) 


9 2-2 in lL> C lias been directly observed through the 
v2 C + 7 -4 3a reaction. 

q Recent ah initio EFT lattice calculations performed by 
Epelbaum, Krebs, Lah.de, Lee, and Meifiner predict the 2.7 
2 MeV above the Hoyle State. 
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Summary 


r 


2 1 



Eres 

E a (res) 

r 70 (res) 

B (E2 : -> G?') 


(MeV) 

(kcV) 

(mcV) 

(e 2 fra 4 ) 

HI7S 

EFT 

10.03(11) 

9.65 

800(130) 

60(10) 

0.73(13) 

2(1) 


9 2t in V2 C has been directly observed through the 
l2 C 4- 7 — > 3a reaction. 

« Recent ah initio EFT lattice calculations performed by 
Epelbanm, Krebs, Lahde, Lee, and Meifiner predict the 2.7 
2 MeV above the Hoyle State. 

■u 

o Revised triple-a reaction rates can affect nucleosynthesis of 
heavy elements during explosive astrophysics scenarios. 
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Success of the two step model at the 2.64 M eV 
3 ■ resonance (r= 400 keV). 
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Cross Section for the HB(p,aJ8Be->aa Reaction 

Although 3 a-particles are emitted, the number in a given energy interval 

depends upon the reaction dynamics. 

Therefore, instead of cross section, we report X = Counts/Luminosity 

ie. X = Counts cm2/sr 
NtNp dW 

Nt = #of target nuclei/cm2 
Np =# of incident protons 
dW = detector solid angle 

X has the same units as differential cross section, but the # of a-particles 
has not been divided out It allows the user to find the number of 
outgoing alpha particles an any energy and angle where X is known for 

any specified Luminosity. 
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